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Abstract 

Background and Objective: Genetic factors are important in the pathogenesis of fractures. Notably, estrogen receptor a 
(ESR1) has been suggested as a possible candidate gene for hip fractures; however, published studies of ESR1 gene 
polymorphisms have been hampered by small sample sizes and inconclusive or ambiguous results. The aim of this meta- 
analysis is to investigate the associations between two novel common ESR1 polymorphisms (intron 1 polymorphisms Pvull- 
rs2234693: OT and Xbal-rs9340799: A>G) and hip fracture. 

Methods: Crude odds ratios (ORs) with 95% confidence intervals (CIs) were used to evaluate the strength of the association. 

Results: Five case-control and three cohort studies were assessed, including a total of 1,838 hip fracture cases and 14,972 
healthy controls. This meta-analysis revealed that the Pvull T allele is a highly significant risk factor for hip fracture 
susceptibility, with an effect magnitude similar in male and pre-menopausal and post-menopausal female patients. In 
stratified analysis based on ethnicity, the Pvull T allele remained significantly correlated with increased risk of hip fracture in 
Caucasian populations; this correlation, however, was not found in Asian populations. Unlike the Pvull polymorphism, we 
did not find significant differences in the Xbal (A>G) polymorphism allele or genotype distributions of hip fracture patients 
and controls. We also found no obvious association between the Xbal polymorphism and hip fracture in any of the racial or 
gender subgroups. 

Conclusion: Our findings show that the ESR1 Pvull T allele may increase the risk of hip fracture and that the Xbal 
polymorphism is not associated with hip fracture. 
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Introduction 

Hip fracture, the most serious complication of osteoporosis, is a 
common musculoskeletal disorder [1]. It accounts for the majority 
of fracture related health care expenditure and mortality in men 
and women over the age of 50 years worldwide [2,3]. The global 
annual rate of hip fracture has been exponentially increasing, with 
a current rate of approximately 1.5 million cases per year [4,5]. In 
addition, with growing aging populations, the number of 
osteoporosis cases is expected to increase to 2.6 million by 2025 
and to 4.5 million by 2050 [6]. Recently, numerous studies have 
attempted to explore the pathogenesis of this disease [7—10]. Bone 
mineral density (BMD) has been found to be an important clinical 
predictor of fracture risk. Most variance in BMD could be due to 
genetic factors, with as much as 65 — 92% of the difference in 
BMD attributable to genetic influences [11,12]. Rapid progress 
has been made in recent years to identify the genes and alleles that 
affect hip fracture risk, such as the vitamin D receptor [VDR) gene, 
insulin-like growth factor I (IGF-I) gene, collagen type I alpha 
\{C0L1A1) gene and estrogen receptor ot gene (ESR1), etc 
[13-15]. 



Estrogen receptor a, a member of the nuclear receptor super- 
family of ligand-activated transcription factors, is one of the key 
mediators of hormonal response in estrogen-sensitive tissues 
[16,17]. The estrogen-ESRl complex is primarily responsible for 
regulating cellular signal pathways in vivo, as well as bone mass in 
skeletal systems [18,19]. It has been shown that serum estradiol 
level may be a predictor of subsequent bone mass density [20,21] 
and risk for osteoporotic fractures [22,23]. Based on the 
observation that osteoporotic fractures are prevalent among 
women, and especially among postmenopausal women, estrogen 
is also believed to play an important role in hip fractures [24]. 
Thus, genetic variations in ESR1, which are caused by alternative 
splicing that alters the expression of ESR 1 , are likely to affect hip 
fracture susceptibility [25]. Several genetic polymorphisms of the 
ESR1 gene, including ESR1 Xbal (rs9340799, A/G) and Pvull 
(rs2234693, C/T), have been investigated for a possible association 
with hip fracture risk [26 — 28]. Two recent genome wide 
association studies (GWAS) have found that several loci in the 
6q25 ESRl region are associated with bone mineral density of the 
hip and spine, providing rapid insights into the genetic association 



PLOS ONE | www.plosone.org 



1 



December 2013 | Volume 8 | Issue 12 | e82806 



ESR1 Gene and Hip Fracture Risk 



between ESR1 gene and hip fracture risk [29,30]. A previous 
meta-analysis by Ioannidis et al. found that the adjusted odds of 
vertebral fracture were reduced by 35% in women homozygous 
for the absence of an Xbal recognition site, while no significant 
effects on fracture risk were observed for PvuII polymorphism 
[31]. However, the results from another recent meta-analysis by 
Lei et al. indicated that a modest, but statistically significant, 
association between the ESR1 PvuII pp genotype and vertebral 
fracture existed in five case-control studies, but no association 
between the PvuII polymorphism and hip fracture were observed 
in one study [32]. 

Two previous meta-analyses have concluded that the ESR1 
Xbal polymorphism may be associated with decreased risk of 
vertebral fracture in women and that the PvuII polymorphism is 
not associated with vertebral fracture risk [31,32]; however, they 
did not provide evidence of the two polymorphisms' correlation 
with hip fracture susceptibility. Gender differences are extensively 
reported in the epidemiology of hip fracture, which is particularly 
prevalent among postmenopausal women. Therefore, we per- 
formed a meta-analysis of all eligible case-control studies 
associating hip fracture risk linked with gender and menopausal 
status in order to explicate the relationship between two common 
polymorphisms (PvuII and Xbal) in the ESR1 gene and hip 
fracture susceptibility. We hope this meta-analysis will help with 
the early identification and therapeutic treatment of hip fracture. 

Materials and Methods 

Literature search 

We performed an electronic search for papers published before 
June 19th, 2013 in PubMed, Embase, Web of Science, Chinese 
National Knowledge Infrastructure (CNKI), and Chinese Bio- 
medical Literature Database (CBM). Literature searches were 
performed by an expert using the following key words and MeSH 
terms: ('fractures, hip' or 'subtrochanteric fractures' or 'femoral 
neck fractures' or 'femur neck fractures'), ('genetic polymorphism' 
or 'single nucleotide polymorphisms' or 'gene mutation' or 'genetic 
variants') and ('estradiol receptor alpha' or 'ER alpha' or 'estrogen 
receptor 1' or 'ESRP). The reference lists of the prospective 
articles were also reviewed to identify additional relevant 
publications and studies. 

Inclusion and exclusion criteria 

Studies included in our meta-analysis had to meet the following 
criteria: (a) population-based case-control, cohort or cross-section- 
al study focusing on associations of ESR1 PvuII and/or Xbal 
polymorphisms with hip fracture risk; (b) hip fracture cases must 
be diagnosed as pathologic fractures or fractures resulting from 
trauma other than a fall; (c) inclusion of sufficient data on sample 
size, odds ratio (OR), and 95% confidence interval (CI); and (d) 
published in the English or Chinese language. 

Studies were excluded when they represented duplicates of 
previous publications, or were meta-analyses, letters, reviews or 
editorial articles. Additionally, when data were included in 
multiple studies using the same case series, either the study with 
the largest sample size or the study with the most recendy 
publication date was selected. All disagreements on study inclusion 
were resolved through discussions. To ensure the rigor of the 
current meta-analysis, we designed and reported the meta-analysis 
according to the Preferred Reporting Items for Systematic 
Reviews and Meta-analyses (PRISMA) statement. The relevant 
checklist is shown in Supplement SI. 



Data extraction 

All data from included studies were extracted independently by 
two investigators using a piloted data standardized form (any 
discrepancies were resolved through discussion and, when 
necessary, adjudicated by a third reviewer): the first author's 
surname, year of publication, country of origin, published 
language, gender and ethnicity of study subjects, study design, 
number of subjects, SNP genotyping methods, genotyping method 
and detected sample, allele and genotype frequencies, and 
evidence of Hardy- Weinberg equilibrium (HWE) in controls. In 
addition, we also compared key study characteristics, such as 
location, publication date and authorship, to determine the 
existence of multiple publications from the same study. 

Quality assessment of included studies 

Two authors independendy assessed the quality of the published 
articles according to the modified STROBE quality score systems 
[33]. Forty assessment items matching with the quality appraisals 
were used in this meta-analysis, with scores ranging from 0 to 40. 
Scores of 0—20, 20—30 and 30 — 40 were defined as low, moderate 
and high quality, respectively. Differences were resolved through 
discussions between the two authors; if no agreement could be 
reached, a third reviewer was consulted. The modified STROBE 
quality score system is available in Supplement S2. 

Statistical analysis 

Crude ORs together with their corresponding 95% CIs were 
used to calculate and assess the strength of associations between 
ESR1 PvuII and/or Xbal polymorphisms and hip fracture risk 
under five genetic models: allele, dominant, recessive, homozy- 
gous, and heterozygous models. The deviation of frequency from 
those expected under Hardy-Weinberg equilibrium (HWE) was 
assessed by Chi-squared goodness of fit tests in the controls. We 
explored inter-study variation by prespecified subgrouping of 
studies according to ethnicity (Caucasian, Asian), gender (female, 
male) and menopausal status among women (pre-menopausal, 
post-menopausal), where applicable. The statistical significance of 
the pooled OR was assessed with a Z test. Between-study variation 
and heterogeneity were estimated using Cochran's (^statistic, with 
P<0.05 as a cutoff for statistically significant heterogeneity [34]. 

We also quantified the effect of heterogeneity with the f test 
(ranges from 0 to 100%), which represents the proportion of inter- 
study variability that can be attributed to heterogeneity rather 
than to chance [35]. The fixed effects model (Mantel-Haenszel 
method) was used, except when a significant Q;test (P<0.05) or 
7 2 >50% indicated the existence of heterogeneity among studies; 
otherwise, the random effects model (DerSimonian- Laird method) 
was applied to the meta-analysis. In order to ensure the reliability 
of results, a sensitivity analysis was performed by omitting 
individual studies. Begger's funnel plots were used to detect 
publication bias. In addition, Egger's linear regression test, which 
measures funnel plot asymmetry via a natural logarithm scale of 
OR, was also used to evaluate publication bias [36]. All P-values 
were two-sided. Analyses were conducted with STATA Version 
12.0 software (Stata Corp, College Station, TX). Meta analysis, 
sensitivity analysis and publication bias were conducted using the 
following STATA's user-written functions: metan, metaninf and 
metabias, respectively. 

Results 

The characteristics of included studies 

Literature search yielded 129 reports, of which five population- 
based case-control [25,26,37-39] and three cohort [27,40,41] 
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studies met the inclusion criteria for studies on the associations 
between hip fracture susceptibility and the two common 
polymorphisms of ESR1, PvuII (rs2234693) and Xbal 
(rs9340799) in intron 1. A flow diagram of the studies selection 
process, as well as the specific reasons for exclusion from the meta- 
analysis is shown in Figure 1 . The selection process involved three 
steps: identification, screening and inclusion. 141 relevant reports 
were initially indentified, 82 were excluded after title and abstract 
review and 51 more were excluded after full-text review. Finally, 
eight studies were included in this meta-analysis. The publication 
years of the selected studies ranged from 2000 to 2012. There were 
six studies on subjects of Caucasian descent and two studies on 
subjects of Asian descent. All included studies extracted DNA from 
peripheral blood and a classic PGR-RFLP assay was used in 6 out 
of 8 studies. Genotype frequencies of healthy controls in all studies 
were consistent with the HWE test. The qualities of the included 
studies were moderately high with a STROBE score greater than 
20. The detailed characteristics of the involved studies are listed in 
Table 1. The genotype distribution and allele frequency data is 
shown in Supplement S3. 

Association between ESR1 PvuII (OT) polymorphism 
and hip fracture risk 

The association of the ESR1 PvuII (OT) polymorphism and 
hip fracture was investigated in eight studies, with a total of 1,838 
cases and 14,972 healthy controls. Among these studies, six were 
conducted in Caucasian populations and two in Asian populations. 



The findings of this meta-analysis on the correlation between the 
polymorphism and hip fracture risk are summarized in Table 2. 

The results revealed that a highly significant increased effect was 
conferred by the PvuII T allele on hip fracture risk, with an 
approximately 18% increment in the odds (for allele model: 
OR=1.19, 95%CI: 1.10-1.28, P<0.001; dominant model: 
OR= 1.25, 95%CI: 1.09-1.44, P=0.001; recessive model: 
OR= 1.26, 95%CI: 1.12-1.42, P<0.001; homozygous model: 
OR= 1.40, 95%CI: 1.19-1.64, P<0.001; heterozygous model: 
OR= 1.21, 95%CI: 1.07-1.36, P= 0.003), and the magnitude of 
the effect was similar in female (for allele model: OR =1.17, 
95%CI: 1.07-1.29, P<0.001; dominant model: OR=1.22, 
95%CI: 1.05-1.42, P= 0.001; recessive model: OR=1.23, 
95%CI: 1.08-1.39, P=0.002; homozygous model: OR=1.35, 
95%CI: 1.14-1.61, P=0.001; heterozygous model: OR=1.18, 
95%CI: 1.03-1.35, P= 0.017) and male (for allele model: 
OR= 1.46, 95%CI: 1.15-1.83, P=0.002; dominant model: 
OR= 1.80, 95%CI: 1.10-2.94, P=0.019; recessive model: 
OR= 1.59, 95%CI: 1.16-2.20, P= 0.005; homozygous model: 
OR = 2.16, 95%CI: 1.27-3.67, P=0.004; heterozygous model: 
OR= 1.44, 95%CI: 1.03-2.20, P= 0.033) (Figure 2B). In 
addition, we isolated the study including pre-menopausal and 
post-menopausal women and analyzed the data in both groups. In 
premenopausal women, we observed that PvuII (C/T) is 
significantly correlated with increased risk of hip fracture under 
allele and homozygous models (OR= 1.30, 95%CI: 1.02-1.66, 
P=0.033; OR= 1.73, 95%CI: 1.04-2.88, P=0.033). In 
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Figure 1. Flow diagram of selection of studies and specific reasons for exclusion from the present meta-analysis. 
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Figure 2. Forest plot of ORs for the association between Pvull (OT) polymorphism and susceptibility to hip fracture in subgroup 
analysis based on ethnicity (A), gender (B), and menopausal status (C) under the allele model. 
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postmenopausal women, we also observed that TT genotype is 
associated with higher risk of hip fracture than either TC or CC 
subjects (for allele model: OR= 1.15, 95%CI: 1.04-1.26, 
P=0.005; recessive model: OR= 1.21, 95% CI: 1.05-1.41, 
P=0.010; homozygous model: OR= 1.31, 95%GI: 1.08-1.59, 
P=0.007; heterozygous model: OR= 1.18, 95%CI: 1.00-1.37, 
P= 0.043) (Figure 2C). Interestingly, in the stratified analysis by 
ethnicity, the PvuII T allele is significantly correlated with 
increased risk of hip fracture among Caucasian populations (for 
aUele model: OR= 1.18, 95%CI: 1.09-1.28, P<0.001; dominant 



model 
model 
model 
model 



OR = 1.26, 95% CI: 
OR = 1.24, 95%CI: 1 
OR= 1.39, 95% CI: 1 
OR= 1.18, 95%CI: 1 



1.09-1.45, P = 0.002; recessive 
10-1.40, P= 0.001; homozygous 
18-1.64, P<0.001; heterozygous 
.04-1.34, P= 0.013), while such 



result was not found in Asian populations (all P>0.05) (Figure 2A). 

Association between ESR1 Xbal (A>G) polymorphism 
and hip fracture risk 

The association between the ESR1 Xbal (A>G) polymorphism 
and hip fracture was investigated only in four studies, with a total 
of 853 cases and 4,522 healthy controls. Among these studies, 
three were on Caucasians and one on Asians. A summary of the 
meta-analysis findings on the association of the polymorphism with 
hip fracture risk is shown in Table 3. The initial meta-analysis 
showed that there was no statistically significant association 
between the Xbal polymorphism and hip fracture under any of 
the genetic models (for allele model: OR = 1.09, 95%CI 



0.97 
0.89 
0.96 
0.93 
0.94 



.23, 
.47, 
.30, 
.57, 
.29, 



P= 0.125; dominant model: OR = 
P= 0.290; recessive model: OR = 
P— 0.161; homozygous model: OR = 
P= 0.167; heterozygous model: OR : 



.14, 
.12, 



95% CI 
95% CI 
1.20, 95%CI 
1.10, 95%CI 



P— 0.245). The association between either ethnicity or 
menopausal status and the polymorphic genotypes were also not 
statistically significant for any of these analyses (all P>0.05, data 
not shown) (Figure 3). 

Sensitivity analysis and publication bias 

Sensitivity analysis was performed to assess the influence of each 
study on the pooled ORs by omitting individual studies. The 
analysis results suggested that no individual study significandy 
altered the pooled ORs for either ESR1 PvuII (C>T) or Xbal 
(A>G) under the allele model (Figure 4), indicating that our 
studies were statistically accurate. 

Begger's funnel plot and Egger's linear regression test were 
performed on the metadata to assess the publication bias of 
individual studies. The shapes of the funnel plots did not reveal 
any evidence of obvious asymmetry for either the ESR1 PvuII 
(OT) or Xbal (A>G) (Figure 5). Egger's test also displayed no 
significant statistical evidence of publication bias (PvuII: t = —0.78, 
P= 0.456; Xbal: t = 0.96, P= 0.440). 

Discussion 

Over the past decade, many linkage and association studies 
have been performed to find candidate genes in order to analyze 
the relationship between genetic factors and complex traits, such 
as BMD and/or fracture risk [42]. It was reported that small 
alterations in BMD could result in a significant difference in 
fracture risk [43]. Thus, many genes in each category were 
proposed for their association with normal BMD variation, 
yielding an ever expanding candidate gene list. One of the most 
widely studied is the ESR1 gene [7]. A genome wide association 
study (GWAS) in 2008 reported that five SNPs (rs9479055, 
rs4870044, rsl038304, rs6929137, rsl999805) in the 6q25 ESR1 
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Figure 3. Forest plot of ORs for the association of Xbal (A>G) polymorphism and susceptibility to hip fracture in subgroup analysis 
based on ethnicity (A) and gender (B) under the allele model. 
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region showed an association with bone mineral density of both 
the hip and spine, suggesting the possible role of the ESR1 gene in 
the pathogenesis of osteoporosis [29]. A more recent GWAS 
further confirmed the association between the ESR1 gene and 
osteoporotic fractures [30]. The most studied variants of the ESR1 
gene are the PvuII and Xbal polymorphisms, which have been 
linked to a lower sensitivity to estrogen [44—46]. However, the 
association between ESR1 gene polymorphisms and fracture risk is 
still controversial and ambiguous. To date, there are only two 
published meta-analyses that evaluate the possibility of a 
significant association between the ESR1 PvuII polymorphism 
and risk of fracture. The meta-analysis by Ioannidiset al. showed 
that the PvuII polymorphism was not associated with fracture risk 
[7], whereas the other one by Lei et al. suggested a modest but 
statistically significant association between the ESR1 PvuII pp 
genotype and vertebral fracture [32]. Therefore, we decided to 



perform a meta-analysis of all eligible case-control studies on hip 
fracture risk in order to reveal a more accurate relationship 
between the PvuII and Xbal polymorphisms of the ESR1 gene and 
risk of hip fracture. 

This meta-analysis evaluated the association between two 
common ESR1 polymorphisms and hip fracture. The results of 
this meta-analysis show that the PvuII polymorphism has a 
statistically significant association with hip fracture, especially in 
Caucasian populations but not in Asian populations. The 
discrepancy between ethnicity subgroups may be due to sample 
size since only 188 cases and 248 controls were studied in Asian 
populations. Other than family history and ethnicity, the most 
well-established risk factors for hip fracture are gender and 
menopausal status in females [47,48]. Therefore, we performed a 
stratified analysis based on gender and menopausal status. Our 
data revealed that the PvuII T allele was a risk factor in the 
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Figure 4. Sensitivity analysis of the summary odds ratio coefficients of Pvull (OT) and Xbal (A>G) polymorphisms are illustrated 
under the allele model. Results were computed by omitting each study in turn. The two ends of the dotted lines represent the 95% CI (A: Pvull; B: 
Xbal). 

doi:1 0.1 371 /journal.pone.0082806.g004 



subgroup of premenopausal and postmenopausal women, and 
men. Unlike the Pvull polymorphism, no significant difference 
was found in the Xbal (A>G) polymorphism allele or genotype 
distribution between hip fracture patients and controls. However, 
this result warrants further investigations since only four studies 
with small sample sizes examined the Xbal (A>G) polymorphism. 
We also performed a stratified analysis based on ethnicity and 
gender for the ESR1 Xbal (A>G) polymorphism. We found no 



obvious association between the Xbal polymorphism and hip 
fracture in either racial or gender subgroups. This discrepancy in 
the associations of the two polymorphisms may be due to sample 
size since there are only eight studies on Pvull and four on Xbal. 
Although these variants alone are not clinically useful in the 
prediction of risk to the individual person, their interactions with 
candidate genes identified by the GWAS may play important roles 
in the pathogenesis underlying hip fracture. 
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Figure 5. Begger's funnel plot of publication bias in the selection of studies for Pvull (OT) and Xbal (A>G) polymorphisms are 
illustrated (A: Pvull; B: Xbal). Each point represents a separate study by the indicated association. Log[OR], natural logarithm of OR. Horizontal 
line, mean magnitude of the effect. 
doi:1 0.1 371 /journal.pone.0082806.g005 



In interpreting the results of this meta-analysis, some specific 
issues should be mentioned. First, as with other complex traits, hip 
fracture risk may be modulated by several other genetic markers 
besides ESR1 gene, including polymorphisms of the VDR gene, the 
COL1A1 gene, the IGF-I gene, and several other candidate 
genes. Thus, our meta-analysis emphasizes that elucidating the 



pathogenesis of fracture demands an investigation into the 
association of many gene variants constituting distinct pathophys- 
iological pathways. Second, we should emphasize that this meta- 
analysis does not address the question of whether the ESR1 gene 
shows evidence of genetic linkage to BMD. Third, we identified 
only two studies from Asian populations and obtained no data 
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from African populations, thus these two racial groups demand 
further studies. Last, this meta-analysis was based on unadjusted 
ORs estimates as a result of the lack of available information, 
preventing a more precise evaluation with adjusted ORs by certain 
covariates such as age, BMI, and smoking status, etc; thus, our 
data is only a conservative estimates of the association between 
ESR1 gene and hip fracture. Ideally, our conclusions will be tested 
by future studies. Despite the limitations listed above, our meta- 
analysis has some strength. To the best of our knowledge, this is 
the first meta-analysis on the relationship between the ESR1 gene 
polymorphisms and hip fracture. We explored inter-study 
variations by subgrouping studies according to ethnicity, gender 
and menopausal status for females. Furthermore, although this 
meta-analysis does not accommodate all previously published 
data, that is limited compared to the evidence here generated. 

In summary, this meta-analysis of eight studies indicates that 
allelic variation in the PvuII (OT) polymorphism of ESR1 gene 
may be a risk factor for hip fracture, with a similar effect 
magnitude in both pre-menopausal and post-menopausal females, 
and male. However, unlike with the PvuII polymorphism, no 
significant difference was found in the Xbal (A>G) polymorphism 
allele or genotype distribution between hip fracture patients and 
controls, either in racial or gender subgroups. Thus, our results do 
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